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Introduction 

This  program  began  May  15,  2007  and  ended  January  14,  2011.  We  have  successfully 
accomplished  not  only  the  goals  stated  in  the  proposal  but  also  beyond.  Fourteen  papers  have 
been  published  or  are  in  press  at  the  Astrophysical  Journal  (ApJ),  Journal  of  Geophysical 
Research  (JGR)  and  in  ASP  Conference  Series  Proceedings  and  three  have  been  submitted  to  Ap. 
J.  and  JGR  journal  and  are  under  review.  Overall,  there  are  seventeen  papers  are  a  result  of  this 
grant.  The  highlights  from  this  program  are  summarized  in  the  following: 

1.  Fragmentation  of  the  Length  of  Strong  Magnetic  Shear  of  the  Main  Neutral  Line  (Lss): 

A  Potential  Sufficient  Condition  for  Solar  Eruption 

A  recent  paper  by  Wu  et  al.  (2009)  has  shown  that  the  well-known  non-potential 
magnetic  field  parameters  (Falconer  et  al.  2002)  are  the  necessary  conditions  but  not  sufficient 
for  the  initiation  of  solar  eruptive  events.  In  this  study,  we  have  used  a  three-dimensional  data- 
driven  MHD  active  region  model  (Wu  et  al.  2006)  which  is  called  Active  Region  Evolution 
Model  (AREM)  to  analyze  the  magnetic  field  structures  of  an  active  region  10720.  By  inputs  of 
the  Big  Bear  Solar  Observatory  (BBSO)  digital  vector  magnetograph  (DGVM)  to  our  AREM 
(Wu  et  al.  2006),  we  simulated  the  non-potential  parameters  (i.e.  Lss,  length  of  strong  magnetic 
shear  of  the  main  neutral  line;  magnetic  flux  (®);  net  electric  current  (In),  and  normalized 
measure  of  the  field  twist)  at  the  time  before  and  after  a  CME  on  January  15,  2005  as  shown  in 
Figure  1 .  The  results  indicate  a  general  characteristic  that  all  four  non-potential  magnetic 
parameters  are  increasing  before  the  solar  eruption,  and  decreasing  after  the  solar  eruption. 

Figure  2  shows  the  detailed  evolution  of  the  length  of  strong  magnetic  shear  of  the  main  neutral 
line  (Lss);  we  found  from  these  results  that  the  eruption  does  not  occur  at  the  left  upper  comer, 
even  the  length  of  Lss  is  very  long,  but  at  the  right  upper  comer  when  the  length  of  the  strong 
magnetic  shear  is  fragmented,  a  CME  was  launched.  This  feature  leads  us  to  believe  that  the 
fragmentation  may  be  an  indicator  for  an  emption  because  it  represents  the  inhomogeneity  of 
magnetic  shear  along  the  neutral  line.  Since  each  segment  has  different  twisting  angle,  this 
indicates  a  spaghetti  type  of  magnetic  field  configuration.  The  MHD  instability  could  easily  be 
triggered,  and  lead  to  the  emption.  Thus,  this  feature  could  be  considered  as  a  “sufficient” 
condition  for  the  solar  eruption.  Most  recently,  we  have  requested  Dr.  David  Falconer  to 
examine  some  old  data  from  NASA/MSFC’s  vector  magnetogram  which  also  exhibits  this  same 
feature. 


Y(0.62arcsecs) 


22:10  22:20  22:30  22:40  22:50  23:00  23:10  23:20  23:30  23:40  23:50 

Time  (UT) 


FIGURE  1.  The  simulated  non-potential  parameters  (i.e.  Lss,  length  of  strong  magnetic  shear  of  the  main  neutral 
line  (unit  2  x  104  km);  magnetic  flux  (<D)  (unit  5  x  1021  Mx);  net  electric  current  (IN)  (unit  3  x  1011  A);  and 
normalized  measure  of  the  field  twist  (a)  (unit  1.5  x  1  0"5  km))  at  various  times  for  AR10720,  January  15,  2005.  The 
horizontal  bold  line  indicates  the  impulsive  phase  of  a  flare  (onset-to-maximum  of  the  soft  X-ray  emission),  the 
vertical  solid  line  shows  the  occurrence  of  a  CME  in  the  LASCO/C2  field  of  view,  and  the  vertical  dashed  line 
indicates  the  back  extrapolated  time  of  the  CME  launch.  (Wu  et  al.  2009) 
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FIGURE  2.  The  simulated  evolution  of  the  length  of  strong  magnetic  shear  of  the  main  neutral  line  (Lss)  of 
AR10720  during  the  period  22:33UT  to  23:15UT  of  January  15,  2005. 


2.  Simulation  of  the  Sun’s  Emission 


Our  newly  developed  data-driven  3D  MHD  model  with  radiation  effects  was  used  to 
reproduce  the  Sun’s  emission.  To  carry  out  this  simulation,  we  used  the  SOLIS  measurements  of 
line-of-sight  (LOS)  magnetic  field  and  transverse  velocity  from  GONG  for  CR2900  as  inputs  to 
the  model  in  which  the  radiative  emission  and  coronal  magnetic  field  are  obtained  as  shown  in 
Figure  3(a,b).  From  this  new  robust  simulation,  we  note  that  the  closed  field  loop  structures 
correspond  to  the  high  intensity  of  the  emission  and  the  open  field  region  exhibits  low  emission 
with  high  speed  solar  wind.  To  verify  the  model,  we  have  compared  the  simulated  emission  (i.e. 
synthetic  image)  of  CR2900  (upper  panel)  and  the  observations  of  CR2900  EIT  304A  image 
(lower  panel)  in  Figure  2.  It  shows  good  agreement. 

Figure  3(a)  shows  the  model  output  of  the  radiative  emission  on  the  photosphere  and  the 
3D  magnetic  field  topology  is  depicted  in  Figure  3(b).  From  these  new  robust  simulations,  we 
note  that  the  closed  field  loop-structure  produces  high  emission  and  the  open  field  region 
exhibits  low  emission  with  high  speed  solar  wind. 


FIGURE  3.  Simulated  (a)  sun’s  emission  on  the  disk  and  (b)  coronal  magnetic  field. 


FIGURE  4.  The  synthetic  image  computed  from  the  model  density  and  temperature  in  unit  DN/s/pixel:  top  panel 
with  GONG’s  observed  transverse  velocity.  The  bottom  panel  is  average  EIT  304  image  during  CR2009  where  “a” 
and  “b”  indicate  the  location  of  the  inside  and  boundary  of  a  coronal  hole. 


3.  Evolution  of  the  Heliospheric  Current  Sheet  (HCS)  and  Coronal  Magnetic  Field 
Configuration  for  Solar  Cycle  23. 

The  structure  of  the  solar  magnetic  fields  varies  systematically  over  a  period  of  about 
eleven  years  in  relation  to  the  four  phases  of  the  solar  activity  cycle.  The  four  phases  of  the  solar 
activity  cycle  are  identified  as:  solar  minimum,  the  rising  phase,  solar  maximum  and  the 
declining  phase.  To  investigate  the  solar  cycle  effects  on  HCS  and  coronal  magnetic  field,  we 
have  employed  our  3D  MHD  model  without  radiation  to  simulate  the  evolution  of  HCS  and 
coronal  magnetic  field  configuration  throughout  Solar  Cycle  23.  This  simulation  is  done  by 
inputting  observed  line-of-sight  photospheric  magnetic  field  data  from  the  Wilcox  Solar 
Observatory  (WSO)  at  the  lower  boundary.  The  simulation  result  is  depicted  in  Figure  5.  Note 
that  the  HCS  at  CR2050  (current  minimum)  is  quite  different  from  CR1910  (i.e.  the  solar 
minimum  at  the  beginning  of  Solar  Cycle  23).  This  observation  shows  the  HCS  at  CR2050 
exhibits  a  high  inclination  which  is  not  flat  as  for  CR  1910.  The  significance  of  this  result  is  to 
provide  the  background  solar  wind  to  study  the  propagation  of  solar  disturbance  which  has  a 
major  effect  on  the  solar  disturbance  induced  shock  arrive  time  at  Earth. 


Figure  5.  Evolution  of  the  HCS  topology  and  coronal  magnetic  field  configuration  throughout  Solar  Cycle  23. 
Five  Carrington  Rotations  are  classified  into  one  group  according  to  their  associated  phases  in  the  solar  cycle.  The 
HCSs  are  shown  above  as  the  isosurfaces  of  Br  =  0  and  the  magnetic  field  configurations  are  shown  below, 
respectively,  for  each  Carrington  Rotation;  Field  lines  have  been  assigned  as  black  for  open  field  lines  and  blue  for 
closed  field  lines;  however,  the  same  foot  points  at  the  photosphere  are  used  in  tracing  the  field  lines  for  each 
Carrington  Rotation  (Flu  et  al.  2008). 
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